Free energy versus internal energy potential for heavy quark systems at finite 

temperature 
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Using the QCD sum rule with its operator product expansion reliably estimated from lattice 
calculations for the pressure and energy density of hot QCD matter, we calculate the strength of 
the J/i/; wave function at origin and find that it decreases with temperature when the temperature 
is above the transition temperature. This result is shown to follow exactly that obtained from the 
solution of the Schrodingier equation for a charm and anticharm quark pair using the free energy 
from lattice calculations as the potential and is in sharp contrast to that using the deeper potential 
associated with the internal energy, which shows an enhanced strength of the J/i/j wave function at 
origin. Our result thus has resolved the long-standing question of whether the free energy potential 
or the internal energy potential should be used in analyzing the spectrum of heavy quark systems 
at finite temperature. 
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Ever since the seminal work by Matsui and Satz [Ij, 
suggesting that the J/'0 suppression could be a signa- 
ture for the formation of the quark-gluon plasma (QGP) 
in rclativistic heavy ion collisions [2| , there have been nu- 
merous theoretical as well as experimental studies from 
SPS to LHC on this subject. On the other hand, it 
is quite discomforting to know that even to this date, 
the question whether J/^ dissolves immediately above 
the critical temperature T(. [3|-|5| or remains bounded at 
higher temperatures p-Q is still not settled. The two 
scenarios have led to different phenomenological models 
based on either complete melting or sequential suppres- 
sions of charmonia in QGP for explaining the observed 
suppression of J /ij: production in rclativistic heavy ion 
collisions, resulting thus in different conclusions about 
the phase transition in QCD and the properties of the 
quark-gluon plasma [9| . 

Given the fact that lattice calculations can calculate 
the heavy quark potential with great precision, it might 
seem that the question could be simply answered by just 
solving the Schrodinger equation and finding the ground 
state eigenvalues. However, it is not so clear whether one 
should use the free energy potential or the internal energy 
potential [lO[. If one uses the free energy as the poten- 
tial, the J /ij} then dissolves slightly above Tc whereas it 
remains bounded up to almost 2Tc if the internal energy 
is used. The difference between the two potentials can 
be traced to the question of whether one should subtract 
out the additional gluonic entropy contribution associ- 
ated with adding the heavy quark pair in the QGP jll| . 
Because of this difference, the internal energy potential 
is deeper and also rises more sharply with increasing in- 
terquark distance, resulting in a J /ip ground state wave 
function that is more localized at origin than that from 



the free energy potential. This effect leads to marked 
different behaviors in the temperature dependence of the 
J /ij} wave function at origin. 

The J /i}) wave function at origin is the non-rclativistic 
limit of the overlap of the J /i]j with the charm vector 
current and is thus a well defined gauge invariant and 
physically observable quantity. In this work, using QCD 
sum rules at finite temperature, we show that the tem- 
perature dependence of this overlap follows exactly the 
temperature dependence of the J /ij} wave function at ori- 
gin that is obtained from using the free energy potential 
in the Schrodinger equation for the charm and anticharm 
quark pair. Unlike the uncertainty of the QCD sum rule 
approaches in predicting particular changes in the whole 
spectral density, the overlap can be obtained from the 
integrated part of the spectral density and thus reliably 
estimated. In fact, we will show that the zero tempera- 
ture sum rule reproduces the strength needed to explain 
the dielectron partial width of the J/fA and that the finite 
temperature calculations are equally reliable. 

We start with the correlation function of the charm 



vector current J,, 
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Inserting J /ipi as an intermediate state using the normal- 
ization {J/'ip{p)\J/ip{p)) = {2Tr)^2ujpS^{p — p), we have 
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where 



(J^(0)|J/V^) = i2^mjf^Nce^,iPiO). 



(2) 
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In the above, e^ is the polarization vector and ■0(0) re- 
duces to the wave function of J /ip at origin in the non- 
relativistic limit. Defining the charm quark pair polar- 



via the following dispersion relation after Borel transfor- 
mation; 



ization function n{q) = ^\li;^{q) / [iq^) , we find the J/V' M{M^) = / dse^'/^^'lmnis), 
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contribution to its imaginary part to be 
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This quantity is related to the dielectron decay width of 



JMr 
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Using the empirical value 



i J, , = 5.55 keV for J/ijj in vacuum 12|, we obtain the 
value 



|V'(0)| = 0.211 ± 0.04 GeV' 



3/2 



(6) 



for the overlap of the charm vector current with a free 
J/ip or the wave function of a free J/ /ip at origin. We 
note that in identifying ip{0) as the J/tp non-relativistic 
wave function at origin, Eg . (|5]) is expected to have sev- 
eral non-trivial corrections [l3|. For example, within the 
non-relativistic limit, we could have used 2mc, with nic 
being the charm quark mass, in the denominator, and 
this would amount to about 20% correction and hence 
the uncertainty given in Eq. 1^. Since the relativistic 
corrections are expected to be temperature independent, 
we will study the temperature dependence of the overlap 
of charm vector current with J/i/i at finite temperature 
relative to its vacuum value. 

The residue of the charmonium correlator /o given in 
Eq. ^ can be reliably calculated in the QCD sum rule 



method [IJ, |15| . While masses of heavy quark systems 



have been mostly studied in the moment sum rule, we use 
here the Borel transformed sum rule, which suppresses 
the contribution from the continuum, to calculate /p as 
it is more robust under changes of the continuum [16|. 
The generalization of the Borel sum rule to finite temper- 
ature is well founded as the operator product expansion 
(OPE) is well known up to dimension 4 operators, which 
are needed for the present analysis, by using the lattice 
data on the energy density and pressure of hot QCD mat- 
ter. The method has been used to study the masses of 
J/0 at finite density [17[ as well as at finite tempera- 
ture [l^l- Recently, the method has been combined with 
the maximum entropy method to reconstruct the char- 
monium spectral density at finite temperature [J]. The 
Borel sum rule for /o only involves the strength at the 
pole, which is an integrated quantity, and thus does not 
depend on the details of the spectral density. Hence, we 
concentrate on the temperature dependence of /q. 

Specifically, the OPE for the charm quark pair polar- 
ization function Tl{q) is equated to the spectral density 



after neglecting the thermal factor tanh[s/2r] in the 
spectral density as its correction to charmonium or the 
continuum is negligible up to temperatures of l.lTc [18|, 
llOf . For the spectral density, we assume the following 
form: 

Imn(s) = Imn^/'^(s) + 9{s - so)ImnP"*(s) (8) 

with the first term given by the pole form in Eq. Q. 

The Borel transformed OPE and its temperature de- 
pendence up to dimension 4 operators are: 

7W(A'/2) ^ e-''nA{i^)[l + asiM^)a{iy) + b{iy)<j)b{T) 

+cMMT)] (9) 

with 1/ = Am'^/M'^ being a dimensionlcss scale parameter. 
The Wilson coefficients a(i^), 6(j/), and c(i^) are summa- 
rized in Ref. [20|, where the temperature dependence of 
the scalar and twist-2 gluon condensates (pb and 0c are 
also given. We note that the truncation of the OPE up 
to dimension 4 operators is valid up to temperature of 
T= l.lTcaswell [H. 

Using the OPE side of Eq.®, we can express the 
residue /o as 
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where the J/0 mass is given by 



V/V' 
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MOPE{AP) - X™nt(^2. 5jj) 



(11) 
The continuum part of the correlator tW^""* depends on 

the threshold parameter sq. Following the method used 
in Ref. [21[ for the estimation of mj/^ , we determine sq 
by requiring |'0(O)| = rrij/^fo/ (12Tr) to be least sensitive 
to M^. We have confirmed that the AP dependence of 
|'0(O)| is quite similar to that of rrij/^, but the resultant 
value of Sq differs slightly from the case in which the 
Borel curve for mj/.^ is optimized. We have also esti- 
mated the systematic error due to the M^ dependence 
of |V'(0)| and found it to be at most 0.001 GeV^/^. We 
have further examined the effects of finite width above 
Tc- Introducing the width, however, brings an ambiguity 
on the stabilization of the Borel curve due to possible 
absence of a solution [20[. Nevertheless, the change due 
to the introduction of nonzero width can be understood 
as a result of the strong correlation among the mass, the 
width and the threshold as shown in Ref. [20|. Since a 
nonzero width leads to a larger mass while keeping the 
residue unchanged, the resultant |0(O)| becomes slightly 



larger as shown in Fig. [l]by the Borel curves for |'!/'(0)| 
at T = 1.0 Tc for three different values of the J/ if) width. 
Nevertheless, the difference is small and our prediction 
for |'!/'(0)| is thus robust. 
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FIG. 1: (Color online) Borel curves for |?/;(0)| at T = 1.0 Tc 
for different values of J/ip width and continuum threshold. 

Figure [5] shows the temperature dependence of |'0(O)I 
and mji^ obtained by the Borel stability analysis at each 
temperature up to slightly above T = 1.05 Tc beyond 
which the OPE becomes less reliable [l8[ . Irrespective of 
the J/ip width used in the analysis, both |'0(O)| (filled 
symbols, left vertical axis) and mj/^ (open symbols, 
right vertical axis) are seen to decrease with increasing 
temperature. 
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FIG. 2: (Color online) Temperature dependence of |i/'(0)| 
(filled symbols, left vertical axis) and mass of J/tp (open sym- 
bols, right vertical axis) obtained from the QCD sum rule for 
different values of J/tp width. 

To find out which potential between a charm and 
anticharm quark pair correctly reproduces |'0(O)I ob- 
tained from the QCD sum rule, we solve the following 



Schrodinger equation between a charm and anticharm 
pair: 



2mr 



1 



V'^ + V{r,T) 



iir,T) = Mi^{r,T), (12) 



where ttIc — 1.25 GeV is the bare mass of charm quark 
and ?/'(r, T) is the charmonium wave function at temper- 
ature T. Introducing the potential energy at infinitely 
large distance, Wr = oo,T), the Schrodinger equation is 
modified to [H [gi] 



— + y(r,T) 

IJlr 



4,{r,T) = -e4>iT,T), (13) 



where V{r,T) = V{r,T) — V{r ~ oo,T), and it van- 
ishes at infinitely large distance, and e = 2toc + V{r = 
oo, T) — Af is the binding energy of J/ijj at temperature 
T. For the heavy quark potential, we use either the free 
energy between a heavy quark-a ntiq uark pair that is ex- 



itiquarK pair tnat is ex- 
M S [23 or the more 



traded from lattice calculations 

attractive one based on the internal energy by adding the 

contribution from the entropy density. 
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FIG. 3: (Color online) Temperature dependence of |i/'(0)| ob- 
tained from the free energy (filled squares) and internal energy 
(filled circles) potentials together with that from the QCD 
sum rule (open circles). 

Figure [3] shows the temperature dependence of the 
strength of the J/ip wave function at the origin obtained 
by solving the Schrodinger with the two different poten- 
tials, together with the result from the QCD sum rule in 
the case of vanishing width. With the internal energy as 
the potential, the strength is seen to increase by almost 
a factor of two at slightly above the critical temperature. 
On the other hand, the strength decreases monotonically 
with temperature when the free energy is used as the po- 
tential, strikingly similar to the result from QCD sum 
rules in both its behavior and values. Even allowing for 



±20% uncertainty in the sum rule result, the case with 
internal energy is outside the range of expectations from 
sum rule calculations. Hence we can conclude that to 
correctly reproduce the non-relativistic wave function of 
J /il) at finite temperature, one should use the free energy 
potential in the Schrodinger equation. 
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FIG. 4: (Color online) Temperature dependence of mj/^ ob- 
tained from solving the Schrodinger equation using either the 
free energy (filled squares) or the internal energy (filled cir- 
cles) and from the QCD sum rule analysis (open circles). 

In Fig. m wc plot the mass of J/ip obtained from 
Eq. P^ with different potentials, and compare the re- 
sults to that obtained from the sum rule analysis of 
Eq. PT|) . Again, we find that the result obtained from 
using the free energy is consistent with that from the 
QCD sum rule analysis. 

The result obtained in the present study that the po- 
tential between a heavy quark-antiquark pair in QGP 
is determined by their free energy has important phe- 
nomenological consequences. For J/ip, this leads to a dis- 
sociation temperature that is only slightly above Tc [25| . 
As a result, all charmonium states would essentially melt 
in the quark- gluon plasma. The observed J/ip^s in rel- 
ativistic heavy ion collisions are then either from initial 
hard collisions that occur outside the QGP or from the re- 
combination of charm and anticharm quarks in the QGP 
during hadronization. 

It should be noted, however, that in obtaining our re- 
sult from the Schrodinger equation, we have kept the 
charm quark mass to its zero temperature value and also 
taken the heavy quark potential to be real. These quanti- 
ties taken separately are gauge dependent, and different 
prescriptions can be used to determine their values at 
finite temperature. Since the heavy quark potential is 
expected to acquire a gauge dependent imaginary part 
at finite temperature [2 a . [27 1 , a stronger potential and 



a different charm quark in-medium mass may then be 
needed to reproduce gauge invariant and physically ob- 
servable quantities like |^(0)|, mj/.^ and the J/ijj disso- 
ciation temperature. 
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